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Abstract 

The increase in moment of inertia of SFg in helium nanodroplets is calculated using the quantum 
hydrodynamic approach. This required an extension of the numerical solution to the hydrodynamic 
equation to three explicit dimensions. Based upon an expansion of the density in terms of the 
lowest four Octahedral spherical harmonics, the predicted increase in moment of inertia is 170uA 2 , 
compared to an experimentally determined value of 310(10) uA 2 , i.e., 55% of the observed value. 
The difference is likely in at least part due to lack of convergence with respect to the angular 
expansion, but at present we do not have access to the full densities from which a higher order 
expansion can be determined. The present results contradict those of Kwon et ah, J. Chem. Phys. 
113, 6469 (2000), who predicted that the hydrodynamic theory predicted less than 10% of the 
observed increase in moment of inertia. 
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There is considerable current interest in the spectroscopy of atoms and molecules solvated 
in liquid Helium, particularly in 4 He nanodroplets []]]. These provide microscopic probes 
of the nature of this unique substance, the only physical system whose equilibrium state 
remains a liquid as T — ► K. Further, helium has properties that make it an almost ideal 
'matrix' for the production and characterization of novel chemical species One such 
property is that even highly anisotropic solutes give rotationally resolved spectra, though 
with effective rotational moments of inertia several times that of the isolated molecule in the 
gas phase H, £|. Development of a quantitative predictive theory for the enhanced effective 
moments of inertia will be valuable for future use of the spectroscopically observed rotational 
constants to help determine the chemical carrier of an unassigned spectral feature, as they 
are often used in gas phase spectroscopy. 

It has been widely recognized that enhanced moments of inertia arise from the kinetic 
energy of helium motion that is correlated with the rotation of the solute. In the case of SF 6 - 
HeTv 0, a fixed frame, fixed node diffusion Monte Carlo calculations on small clusters with 
N = 8 — 20 He atoms have recovered rotational excitation spectra in excellent agreement 
with that observed for SFq in He nanodroplet with much larger number of helium atoms. 
Such calculations do not directly give any dynamical information, however, and thus leave 
open the question of how to physically characterize the helium motion. Two quantitatively 
predictive dynamical models have been put forth that invoke very distinct types of helium 
motion H, H. In one [H, [7J, a proposed 'nonsuperfluid fluid' density of helium was calculated 
using Path Integral Monte Carlo (PIMC) methods, and this density is assumed to rotate 
rigidly with the molecule provided that the molecule-He potential is sufficiently anisotropic 
relative to the induced rotational energy. To date, this two fluid model has been applied only 
to SF 6 and OCS || [7j. For SF 6 , the 'two fluid model' gave a result in excellent agreement 
with experiment ||. The predicted increase in the effective moment of inertia, AI, was 
calculated as 327uA 2 compared to a valued of 310(10) uA 2 calculated from the observed 
rotational constants. The level of agreement, fact, is far better than should be expected 
given the considerably uncertainty in the relevant He-SF6 interaction potential 0. The two 
parameters that determine the anisotropy of the He-SF 6 potential well have reported values 
of -0.6 ± 0.3 and 0.14 ± 0.14 in Ref. §. 

The second approach, published by ourselves and coworkers M, is based upon a hydro- 
dynamic treatment for the helium flow, which is assumed to maintain a constant solvation 
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density in the frame rotating with the molecule and be ideal (aviscous and irrotational). The 
hydrodynamic approach was applied to a number of linear molecules, and was found to be in 
good agreement for heavier molecules, including OCS. A key assumption of this work is that 
of adiabatic following by which we mean that the helium density in the molecular frame of 
a rotating molecule is the same as for a rigid molecule, fixed in the laboratory frame. It has 
recently been established experimentally jllj that for lighter molecules, in particular HCN 
and DCN, that adiabatic following breaks down and explains at least the sign of the error 
of the hydrodynamic theory in these cases. 

In order to apply the hydrodynamic model, one needs the three dimensional helium den- 
sity around the solute, which is known to be highly structured. The earlier hydrodynamics 
work used helium density functional theory (DF) [11] to estimate this quantity. For sys- 
tems with cylindrical or higher symmetry, DF is many orders of magnitude computationally 
less expensive than Quantum Monte Carlo Methods but introduces additional uncertain- 
ties beyond the ever present uncertainty in the Helium-Solute interaction potential. In this 
work, the helium density around a static SFq molecule, previously published by Barnett and 
Whaley is used to calculate the hydrodynamic contribution to the moment of inertia of 
this molecule. Because of the different symmetry (Oh versus C^) of this solute compared 
to those treated previously, some changes had to be made in the computational procedure. 
Most significant, one degree of freedom could not be removed by separation of variables 
(as it could be for a cylindrically symmetric density) and, as a result, the hydrodynamic 
equation for the velocity potential had to be solved numerically in three degrees of freedom. 
The computational procedures used in the present case are described in the next section and 
the results presented in the last section. 



I. HYDRODYNAMIC CALCULATIONS 



The hydrodynamic calculations require an estimate of the three dimensional helium den- 
sity around the solute molecule, SF 6 in this case. This density will be totally symmetric in 
the point group of the molecule, Oh in this case. This density, p(r, Q), can be expanded in 
terms of the spherical tensor operators: 

p(r, Q) = p (r) + a 4 (r) T 4 (fi) + a 6 (r) T 6 {Q) + a 8 (r) T 8 (f2) + . . . (1) 
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where Q = (8, if) and Tl are linear combinations of the spherical harmonics, Ylm's, that 
transform as Ai g in the Oh point group. In particular: 

t 4 = y^y 40 + (y 40 + f 4 _ 4 ) (2) 

T 6 = ^F 60 - ^ (Feo + F 6 _ 4 ) (3) 

T 8 = ^Y m + ^ (F 84 + F 8 _ 4 ) + (F 88 + F 8 _ 8 ) (4) 

Fox and Ozier [O] have presented a general procedure for to calculate these and higher 
harmonics. 

Barnett and Whaley have calculated the helium density around SF 6 using the diffusion 
Monte Carlo method. Their paper contains figures which give the radial isotropic density, 
Po(r), and radial 'cuts' of the density along the C 2 , C 3 , and C 4 symmetry axes. If we assume 
that the density contains only the terms explicitly given in Eq. [I], then we can use the 
isotropic density and the density cuts to determine the radial coefficients a^r) , a e (r) , and 
a 8 (r) by using the equations: 

(5) 
(6) 
(7) 
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In these equations, pi represents the radial density along the Q symmetry axis. The radial 
density functions where determined by digitization of the images from the paper by Barnett 
and Whaley for clusters of 69 Helium atoms, the largest for which they reported the density 
cuts. These equations were derived by evaluation of the tensor operators (Tj) along each of 
the symmetry axes and then inverting the linear system that relates the calculated density 
along each cut. Figure |l] shows the radial tensor densities calculated from these density cuts. 
It is apparent from this figure that the density in the first solvation shell (which is the most 
important for determination of the moment of inertia) is likely not fully converged with the 
present truncation of the density expansion. Unfortunately, we were not able to obtain the 
primary data from reference which would have allowed the expansion to be carried to 
higher order. 
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The assumption that the helium motion is irrotational (V x v = 0) implies that the 
velocity, v, can be written as the gradient of a scalar function, v = — V0, where is known 
as the velocity potential. In order to calculate the increased moment of inertia caused by 
motion of the helium, we will assume that the molecule is undergoing classical rotation 
with angular velocity uo around the z axis. The equation of continuity gives the following 
equation: 



where the second equality comes from the assumption that the helium density remains 
constant in the frame rotating with angular velocity to around the z axis. The boundary 
conditions on are that the normal component of the fluid velocity match the normal 
component of velocity of any surface. The solution to Eq. [8] is linear in the angular velocity 
and thus is solved numerically for unit velocity. Equating the helium kinetic energy with 
the rotational kinetic energy, we can define A/h by: 



The first equation holds for any velocity potential. The second has been derived from the 
first using vector identities and Eq. ||. As such, these two estimates need only be equal for 
that is a solution to Eq. |8|. It can also be shown that the net orbital angular momentum 
produced by the helium flow is, for the solution of Eq. ||, Jn e = AI h u}. 

A numerical solution for was found on a grid of points in a spherical coordinate system. 
The three cartesian axes of this coordinate system is aligned with three of the S-F bonds. 
The range of r was selected between 3.5 and 10 A with between 51 and 201 radial points. 
The normal component of was selected to be zero on both the inner and outer constant 
r surfaces. must be invariant to reflection in the x, y plane, and thus we numerically 
restricted the solution to the 9 = [0, 7r/2] and restrict the solution to have zero normal 
derivative on the x, y plane. Reflection in each of the four planes perpendicular to the z axis 
is equivalent to change in the direction of rotation, and thus leads to a change in the sign of 
0. As a result, the numerical solution for could be restricted to the domain tp = [0, 7r/4] 
and selected to be zero on each of the planes ip = 0and7r/4. Between 41 and 161 equally 
spaced angular points were used in each of the angular coordinates coordinates. Since or 
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its normal derivative is equal to zero on all the boundary surfaces, the surface integral in 
Eq. is zero. 

The inhomogeneous partial differential Eq. |8] was converted to a finite difference equation 
in the grid of points and solved by Gauss-Seidel iteration with successive overrelaxation JI1 



These equations (which involve In p and its gradient) become singular when p = 0. As a 
result, the density is bounded to remain above a threshold value, which was selected as 
Pmin = 10~ 5 A~ 3 in this work. Solutions were iterated until the mean squared change in (j) on 
the grid points was less than a fixed fraction (10~ 10 in this work) of the mean squared value 
of (j> on the grid points. For the largest grid used, this required about 6000 iteration cycles 
with an overrelaxation parameter of 1.4 (u in the notation used in Numerical Recipes [14]). 
The two integral estimates for AI^ are found to be 171 and 168 uA 2 respectively. This can be 
compared with a value of AI = 310(10) uA 2 inferred from the observed effective rotational 
constant of SF 6 in 4 He nanodroplets 0. 

Kwon and Whaley |J reported the 'superfluid' densities along the same symmetry axes, 
as calculated using PIMC and their proposed superfluid estimator. These have greatly de- 
creased anisotropy compared to the total density. Repeating the hydrodynamic calculation 
using these densities gave AI^ = 12 uA 2 . This can be compared to their 'superfluid' con- 
tribution of the moment of inertia increase of 22 uA 2 J7|. These authors also claim to have 
calculated the hydrodynamic contribution to the rotational constant used the total density 
around SF 6 and report AI^ = 31 uA 2 . This value is 5.5 times smaller than the value we have 
calculated. Based upon our extensive experience with similar hydrodynamic calculations, 
the modest 50% in hydrodynamic moment of inertia is inconsistent with the considerably 
greater angular anisotropy of the total density compared to their reported 'superfluid' den- 
sity. 



II. DISCUSSION 



The hydrodynamic prediction for the enhanced moment of inertia of SFg in helium is 
~ 55% of the experimentally observed value. This can be contrasted with the results reported 
in [§], where for heavier rotors, the theory appeared to systematically overestimate the 
size of the increased moment of inertia. Similarly, an overestimate of the hydrodynamic 
estimate for the increased moment of inertia was found for a model problem of a planer 
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rotor interacting with a rigid ring of He atoms [[HJ. One explanation for the underestimate 
in the present case is that the density model used has underestimated the true anisotropy of 
the helium density, particularly in the first solvation shell where the hydrodynamic kinetic 
energy density is highest. This is supported by the results presented in fig 1, where it is 
demonstrated that the highest order anisotropy we have retained is still of considerable size 
in this first shell. Since the (p derivative of the density is the source term in the hydrodynamic 
equation, addition of higher order anisotropics is expected lead to increased hydrodynamic 
motion and thus an increased estimate of AI^. This explanation for the limited success of 
the hydrodynamic model in the present case can be tested by repeating the hydrodynamic 
calculation using the full anisotropic helium density calculated by Quantum Monte Carlo 
methods. 

An alternative explanation, which cannot be ruled out at present, is that some fraction 
of the helium motion is not 'irrotational'. The Kelvin minimum energy principle |16| states 
that any 'rotational' solution to the equation of continuity will generate higher helium ki- 
netic energy and thus He contribution to the moment of inertia. This would include the 
hydrodynamics assumed by the 'two fluid' model of Kwon and Whaley |J. 

The fixed node, Diffusion Monte Carlo calculations of Lee et were in excellent 

agreement with the experimental rotational excitation energies. This suggest that, barring 
cancelation of errors, that the nodal structure of the wavefunction assumed in that work 
(which had the nodal properties of the rigid rotor wavefunction for the SFg) should be a 
reasonable description of the true many body wavefunction for this system. It is useful, 
therefore, to examine how consistent the hydrodynamic model is with the functional form 
assumed in that work. 

In a future publication, we will present a quantum derivation of the hydrodynamic ap- 
proach. In that work, it will be shown that Eq. § arises from a variational optimization of 
a one particle phase function that multiples the ground state wavefunction that describes 
the helium in the frame that rotates with the molecule. The wavefunction for the orien- 
tation of this axis system in space will be the rigid rotor function for the molecule. Thus, 
the presence of the one particle phase functions will modify the nodal surfaces from that 
assumed in the Lee et al. calculation. The size of the one particle phase argument will be 
m-R e (f)/h. The maximum value for </> for our solution is ~ 5A 2 uj, and for the J = 1 level we 
can approximate u = 4ttB c s, where _B e ff is the effective moment of inertia of SFq in helium, 
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1.04 GHz. This gives a maximum hydrodynamic phase of 0.04, and thus it appears that 
this will lead of small changes in the nodal properties. We would like to point out that for 
the model problem of the planer rotor coupled to a ring of Helium (which can of course be 
solved exactly), the errors in the equivalent fixed node approximation was of this same size 
and yet the DMC estimate of the rotational excitation energy had errors of at most a few 
percent. 
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FIG. 1: The tensor densities, po(r), a±{r), ae(r), and ag(r), in A , as a function of the He-SFg 
radial distance in A 
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